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Abstract. We present a survey of the formaldehyde emission in a sample of eight Class protostars obtained with the IRAM 
and JCMT millimeter telescopes. The range of energies of the observed transitions allows us to probe the physical and chemical 
conditions across the protostellar envelopes. The data have been analyzed with three different methods with increasing level of 
sophistication. We first analyze the observed emission in the LTE approximation, and derive rotational temperatures between 
11 and 40 K, and column densities between 1 and 20 x 10 13 cirT 2 . Second, we use a LVG code and derive higher kinetic 
temperatures, between 30 and 90 K, consistent with subthermally populated levels and densities from 1 to 6 x 10 s cirr 3 . 
The column densities from the LVG modeling are within a factor of 10 with respect to those derived in the LTE approximation. 
Finally, we analyze the observations based upon detailed models for the envelopes surrounding the protostars, using temperature 
and density profiles previously derived from continuum observations. We approximate the formaldehyde abundance across the 
envelope with a jump function, the jump occurring when the dust temperature reaches 100 K, the evaporation temperature of 
the grain mantles. The observed formaldehyde emission is well reproduced only if there is a jump of more than two orders of 
magnitude, in four sources. In the remaining four sources the data are consistent with a formaldehyde abundance jump, but the 
evidence is more marginal (< 2 <x). The inferred inner H 2 CO abundance varies between 1 x 10~ 8 and 6 x 10~ 6 . The absolute 
values of the jump in the H2CO abundance are uncertain by about one order of magnitude, because of the uncertainties in 
the density, ortho to para ratio, temperature and velocity profiles of the inner region, as well as the evaporation temperature 
of the ices. We discuss the implications of these jumps for our understanding of the origin and evolution of ices in low mass 
star forming regions. Finally, we give predictions for the submillimeter HiCO lines, which are particularly sensitive to the 
abundance jumps. 
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1. Introduction 

Low mass protostars form from dense fragments of molecular 
clouds. During the pre-collapse and collapse phases, the physi- 
cal and chemical composition of the matter undergoes substan- 
tial, sometimes spectacular, changes. From a chemical point 
of view, the pre-collapse phase is marked by the freezing of 
molecules onto the grain mantles. In the very inner parts of 

the pre-stellar condensations, molecules have been observed to 

■ ■ 1 

progressively disappear from the gas phase (e.g. Tafalla et al. 
l2002HBergin et alJl2002h . The CO molecule, whose condensa- 
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tion temperature is around 20 K, is the best studied species both 
because it is the most abundant molecule after H2, and because 
of its important role in the gas thermal cooling. CO depletion 
of more than a factor of ten has been observed i n the centers of 
these condensations JCaselli et alJll998ll2002t bacmann et alJ 
I2002I) . This large CO depletion is accompanied by a variety 
of changes in the molecular composition; the most spectacu- 
lar is the dramatic increase in the molecular deuteration (up to 
eight orders of magnitude with respect to the D/H elementa l 
abundance) observed in formaldehyde (Bacma nn et alj E003). 
The changes are recorded in the grain mantles, where the pre- 
collapse gas will be progressively stored. When a protostar is 
finally born, the dust cocoon warms up and the mantle species 



2 



S. Maret et al.: The H2CO abundance in low mass protostars 



evaporate into the gas phase, returning information from the 
previous phase. 

Most of the studies of the composition of the grain mantles 
have been so far carried out towards massive protostars, be- 
cause they have strong enough IR cont inua against which the 
absorption of ices can be observed ( e.g. lGerakines et alj fl999: 
iDartois et alJl999HGibb et alfeOOCl) . The absorption technique 
allows one to detect the most important constituents of the 
grain mantles: H?Q, CO, CO?, and sometime N H^, CH3OH and 
H 2 CO JSchutte et alJll996HKeane et alJl200lh . In much cases, 
the mantle composition of low mass protostars has been di- 
rectly observed. In these cases, the observations have been car- 
ried out towards protostars that possess a strong enough IR con- 
tinuum (e.g. lBoogert et alJEo OOb). If our understanding of the 
evolution of a protostar is basically correct, those protostars, 
typically Class I or border line Class II sources, represent a rel- 
atively evolved stage, wher e most of the original envelope ha s 
already been dispersed (e.g. lShu et al.ll987llAndre et alJ20 00). 
Furthermore, the observed absorption may be do minated by 
foreground molecular clouds (Booge rt et alJl2002h . Thus, di- 
rect observations of the chemical composition of the primeval 
dust mantles of low mass protostars have so far proven to be 
elusive. 

Alternatively, one can carry out an "archeological" study, 
looking at the composition of the gas in the regions, which 
are known or suspected to be dominated by the gas desorbed 
from grain mantles. This technique has the advantage of be- 
ing much more sensitive than the absorption technique, as it 
can detect molecules whose abundance (with respect to H2) 
is as low as ~ 1CT 11 against a limit of ~ 1CT 6 - 10~ 7 reach- 
able with the absorption technique. Indeed, several very com- 
plex molecules observed in the warm (> 100 K) gas of the 
so called hot cores have been co nsidered hallmarks of grain 
mantle evaporation products (e.g. Bla ke et al.l ll987). Once in 
the gas phase, molecules like formaldehyde and methanol, ini- 
tially in the grain mantles, trigger the formation of more com- 
plex molecules, referred to as daughter or second-gene ration 
molecules (e.g. ICharnlev et al.lll992l ICaselli et alJll993h . The 
gas temperature and density are other key parameters in the 
chemical evolution of the gas, which has the imprint of the pre- 
collapse phase. 

So far, hot cores have been observed in massive protostars, 
and are believed to represent the earliest stages of massive star 
formation, when the gas is not yet ionized by the new born star 
(Kur tz et al.M 2000). Recently, however, it has been proposed 
that low mass protostars might also harbor such hot cores. Note 
that the definition of hot core is not unanimous in the litera- 
ture. Here we mean a region where the chemical composition 
reflects the evaporation of the ice mantles and subsequent reac - 
tions between those species (e .g. iRodgers & Charnlevll2003l) . 
In this respect, Ceccare fh et alJ ll2000alcl) claimed that the low 
mass protostar IRAS 16293-2422 shows evidence of an in- 
ner region (~ 400 AU in size) warm enough (> 100 K) to 
evaporat e the gr a in ma ntles, a claim substantially confirmed 
bv ISchoier et al j (|2002). Indee d, very recent observ ations by 
Cazaux et al. J2003t see also ICeccarelli et al1l2000bl) reveal 
also the presence of organic acids and nitriles in the core of 
IRAS 16293-2422, substantiating the thesis of a hot core region 



in which not only the ices have evaporated but a lso a subse- 
quent hot core chemistry has ensued. Furthermore. iMaret et alJ 
(2002) argued that NGC1333-IRAS4, another low mass very 
embedded protostar, has also such a warm region, somewhat 
less than 200 AU in size. 

Formaldehyde is a relatively abundant constituent of the 
grain mantles and it is a basi c organic molecule t hat forms 
more complex molecules (e.g. ICharnlev et aT1ll992h . For this 
reason, we studied the formaldehyde line emission originating 
in the envelopes of a sample of very embedded, Class low 
mass protostars. In this article we report the first results of this 
systematic study. This is part of a larger project aimed to char- 
acterize as far as possible the physical and chemical composi- 
tion of low mass proto stars during the first phases of formation. 
Jol'g ensen et al.l(Eo 02) determined the temperature and density 
structure for these sources and the CO abundance in the outer 
regions. A forthcoming paper will address the methanol line 
emission in the same source sample, as methanol is another 
key organic mantle constituent, linked by a common formation 
route with formaldehyde. 

One of the ultimate goals of the present study is to under- 
stand the efficiency of H2CO against CH3OH formation in low 
mass protostars, whether and how it depends on the parental 
cloud, and to compare it with the case of massive protostars. An 
immediate goal of the present article is to study the formalde- 
hyde abundance profile in the surveyed sample of low mass 
protostars. In a previous study that we carried out towards 
IRAS 16293-2422, we concluded that formaldehyde forms on 
grain mantles and is trapped mostly in H20-rich ices in the in- 
nermost regions of the envel ope and mostly in CO-ric h ices in 
the outermost regions JCeccarelli et alJEoOOd l200ll) . As the 
dust gradually warms up going inwards, formaldehyde is re- 
leased from the icy mantles all along the envelope. In the hot 
core like region (r < 200 AU) the formalde hyde abundance 
jumps by about a factor 1 00 to ~ 1 x 10 -7 JCeccarelli et alJ 
2000c; ISchoier et alJl2002l) . Similarly, formaldehyde enhance- 
ment is observed in seve ral outflows, because of ice man- 
tle sputtering in shocks dBachiller & Perez Gutierre^ Il997l 
iTafalla et alJ2000l) . In contrast, no jump of formaldehyde abun- 
dance has been detected in t he sa mple of massive protostars 
studied bv Ivan der Tak et all d2000l) . To firmly assess whether 
and by how much formaldehyde is systematically more abun- 
dant in the interiors of low compared to high mass proto- 
stars, a survey of more low mass protostars has to be carried 
out. This will allow us to answer some basic questions such 
as how, when and how much formaldehyde is formed on the 
grain mantles. Given that it forms more complex molecules 
(e.g. Bernstein et alJ 1999) knowing the exact abundance of 
formaldehyde is fundamental to answer the question of whether 
or not pre- and/or biotic molecules can be formed in the 200 or 
so inner AUs close to the forming star. 

In this article we report observations of formaldehyde emis- 
sion in a sample of eight Class sources. After a preliminary 
analysis (rotational diagrams and LVG analysis), the observa- 
tions are analyzed in terms of an accurate model that accounts 
for the temperature and density gradients in each source, as 
well as the radiative transfer, which includes FIR photon pump- 
ing of the formaldehyde levels. The article is organized as fol- 
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lows: we first explain the criteria that lead to the source and 
line selection and the observations carried out (§2). In §3 we 
describe the results of the observations, in §4 we derive the 
approximate gas temperature, density and H2CO column den- 
sity of each source by means of the standard rotational diagram 
technique and by a non-LTE LVG model. In §5 we derive the 
formaldehyde abundance in the inner and outer parts of the en- 
velope of each source, with an accurate model that takes into 
account the structure of the protostellar envelopes. Finally, in 
§6 we discuss the implications of our findings, and conclude in 
§7. 

2. Observations 

2.1. Target and line selection 

We observe d a sample of eigh t protostars, all of them Class 
sources dAndre et alJ l2000l) located in the Perseus, p- 
Ophiuchus and Taurus complexe s, except LI 157-MM th at lies 
in an isolated clump (Bachiller & Perez Gutierrea ll997). The 
eight selected sources ar e among the brightest Class sources 
in the lAndre et alJ l2000) sample. Their physical structure (dust 
density and temperature profiles ) has been determ ined from 
their continuum emission by Jcir gensen et al.1 J20 02) except for 
L1448-N, which is analyzed in this pap er (see AppendixlAl. 

The source distances quoted by J0rgensen et al 
were adopted. The list of the selected sources is reported 
in Table ^ together with their bolometric luminosity, enve- 
lope mass, the ratio of the submillimeter to bolometric lumi- 
nosity, and the bolometric temperature and distance. In the 
same Table, we also report the data relative to IRAS 16293- 
242 2, which was previously obse rved in H?CO trans itions 
bv Ivan Dishoeck et alJ (1 1995b and Loinard et al.1 d2000t). an d 
studied in ICeccarelli et all OOOOd) and lSchoieretal.H 2002). 
IRAS 16293-2422 will be compared to the other sources of the 
sample. 

A list of eight transitions was selected, three ortho-H2CO 
transitions and five para-H2CO transitions (Table |2}. When 
possible, the corresponding isotopic lines were observed in or- 
der to derive the line opacity. The transitions were selected 
to cover a large range of upper level energies (between ~ 20 
and ~ 250 K) with relatively large spontaneous emission co- 
efficients A u j (> 10 -4 s _1 ). The latter condition is dictated 
by the necessity to detect the line, whereas the first condi- 
tion aims to obtain lines that probe different regions of the en- 
velope. Finally, practical considerations, namely having more 
than one line in a single detector setting, provided a further con- 
straint. In the final sel ection we we re helped by our pilot study 
on IRAS 16293-2422 dCeccarelli et alJ EoOOc) and by previous 
studies of the formald ehyde emission in molecular clouds and 
protostellar envelop es (Mang um & Woottenll993Hjansen et all 
ll994ll995l:rCeccarelli et al.l2003l) . 

The formaldehyde transitions between 140 and 280 GHz 
were observed with the single dish IRAM-30m telescope 1 , lo- 
cated at the summit of Pico Veleta in Spain. Higher frequency 

1 IRAM is an international venture supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 
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lines were observed at the JCMT 2 , a 15 meter single dish tele- 
scope located at the summit of Mauna Kea, Hawaii. The choice 
of the two telescopes allows us to have roughly similar beam 
sizes over the observed frequencies. 

2.2. IRAM observations 

The IRAM observations were carried out in November 1999 3 
and September 2002. The various receivers available at the 
30-m were used in different combinations to observe at least 
four transitions simultaneously. The image sideband rejection 
was always higher than 10 dB, and typical system tempera- 
tures were 200-300 K at 2 mm, and 200-500 K at 1 mm. The 
intensities reported in this paper are expressed in main beam 
temperature units, given by: 

r mb = (i) 

where B e g is the main beam efficiency, and F e ff is the forward 
efficiency. The main beam efficiency is 69%, 57% and 42% at 
140, 220 and 280 GHz respectively, and the forward efficiency 
is 93%, 91% and 88% at the same frequencies. Each receiver 
was connected to an autocorrelator unit. For the 1 and 2 mm 
bands, a spectral resolution of 80 kHz and a bandwidth of 80 
MHz was used. These spectral resolutions correspond to a ve- 
locity resolution of 0.09-0.17 km s _1 depending on frequency. 
All IRAM observations were obtained in position switching 
mode. The absolute calibration was regularly checked and was 
about 10%, 15% and 20% at 140-170 GHz, 220-240 GHz and 
280 GHz respectively. Pointing was also regularly checked and 
was better than 3". 

2.3. JCMT observations 

The JCMT observations were obtained in February 2001, 
August 2001 and February 2002. The single sideband 
dual polarization receiver B3 was used with the Digital 
Autocorrelation Spectrometer (DAS). Typical system temper- 
atures were 400 to 800 K. A spectral resolution of 95 kHz for 
a 125 MHz bandwidth was used for most of the lines, while a 
resolution of 378 kHz for a bandwidth of 500 MHz was used to 
observe some of the lines simultaneously. These spectral res- 
olutions correspond to a velocity resolution of 0.08-0.32 km 
s _1 . The antenna temperatures were converted into main beam 
temperature scale using 4 

T mb = -± (2) 

?7mb 

2 The JCMT is operated by the Joint Astronomy Center in 
Hilo, Hawaii on behalf of the present organizations: The Particle 
Physics and Astronomy Research Council in the United Kingdom, 
the National Research Council of Canada and the Netherlands 
Organization for Scientific Research. 

3 IRAM N ovember 1999 data have also been presented in 
lLoinardetal]j2002iJ) 

4 JCMT does not follow the same convention for antenna temper- 
ature than IRAM, the JCMT antenna temperature being already 
corrected for the forward efficiency of the antenna. 
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Table 1. The observed sample. IRAS 16293-2422, which has been studied elsewhere (see text) is included for comparison. 



Source 


a(2000) 


5(2000) 


Cloud 


Dist. 11 


I b 


M env b 


T I T C 

■^smiTi / "^bol 


T b 
1 bol 










(pc) 


(Lo ) 


(Mo ) 


(%) 


(K) 


NGC1333-IRAS4A 


03:29:10.3 


+31:13:31 


Perseus 


220 


6 


2.3 


5 


34 


NGC1333-IRAS4B 


03:29:12.0 


+31:13:09 


Perseus 


220 


6 


2.0 


3 


36 


NGC1333-IRAS2 


03:28:55.4 


+31:14:35 


Perseus 


220 


16 


1.7 


< 1 


50 


L1448-MM 


03:25:38.8 


+30:44:05 


Perseus 


220 


5 


0.9 


2 


60 


L1448-N 


03:25:36.3 


+30:45:15 


Perseus 


220 


6 


3.5 


3 


55 


L1157-MM 


20:39:06.2 


+68:02:22 


Isolated 


325 


11 


1.6 


5 


60 


LI 527 


04:39:53.9 


+26:03:10 


Taurus 


140 


2 


0.9 


0.7 


60 


VLA1623 


16:26:26.4 


-24:24:30 


p-Ophiuchus 


160 


1 


0.2 


10 


<35 


IRAS16293-2422 d 


16:32:22.7 


-24:38:32 


p-Ophiuchus 


160 


27 


5.4 


2 


43 



a Froml Andre et aljf2 000\ except for Perseus sources JCemisll99Cl) . 

From l0raensen et al . (2002) 
From lAndre et alJfeOOd) 
d Included for comparison. 



where rj^ is the main beam efficiency, equal to 63% at the ob- 
served frequencies. The calibration and pointing were regularly 
checked using planets and were found to be better than 30% 
and 3" respectively. The JCMT observations were obtained in 
beam switching mode, with a 180" offset. 



3. Results 

The observed H2CO line spectra are shown in Figs.Hand|2]and 
the results of the observations are summarized in Table|2] 

Most of the lines are relatively narrow (FWHM ~ 2-3 km 
s~') with a small contribution (< 5%) from wings extending 
to larger velocities. The higher the upper level energy of the 
transition the lower the contribution of the wings, which prac- 
tically disappear in the lines observed with JCMT. NGC1333- 
IRAS4A and NGC1333-IRAS4B represent an exception to this 
picture. The line spectra of these two sources are broad (~ 5 
km s~') and the wings are more pronounced than in the other 
sources. Evidence of self-absorption and/or absorption from 
foreground material is seen in most sources, in particular in 
low lying lines. 

In this study we focus on the emission from the envelopes 
surrounding the protostars. Hence, we are interested in the in- 
tensity of the narrow component of the lines, that we fitted 
with a Gaussian. In some cases, a residual due to the "high" 
velocity wings remains, and that has not been included in the 
line flux estimate. When a Gaussian fitting was not possible 
because of self-absorption, the flux of the lines was estimated 
by integrating over a velocity range of + 2 km s around the 
source V] sr . This velocity range corresponds to the width of the 
lines with high upper level energies, where self-absorption is 
less important. For these lines, the self-absorption is included 
in the line flux determination, and the flux measured is there- 
fore slightly smaller than the one that would have been obtained 
by a Gaussian fitting. We observed the brightest lines in the 13 C 
isotopomer of formaldehyde, as reported in Table [5] 



Finally, the errors quoted in Table|2]and[3]include both the 
statistical uncertainties and the calibration error. For non de- 
tected lines we give the 2cr upper limit defined as follows: 



Fmax = 2(1 +a)cr VAv 6v 



(3) 



where <x is the RMS per channel, Av is the line width estimated 
from detected lines on the same source, 5v is the channel width, 
and a is the calibration uncertainty. 

4. Approximate analysis 

4.1. Line opacities 

The detection of some H^CO transitions allows to estimate 
the line opacities of the relevant HpCO lines. Using the escape 
probability formalism and assuming that the H^CO lines are 
optically thin, the ratio between the H^CO and H^CO line 
fluxes can be expressed as: 



"h' 2 co 

7 H»C0 



[ 12 C] 
[ 13 C] 



78 



(4) 



where /3 is the escap e probability, which, in the case of a homo- 
geneous slab of gas Jde Jong et alJl9801) . is equal to: 



1 - exp(-3r) 



3t 



(5) 



and 



jr^j is the is otopic elemental ratio, equal to 70 
(Boo gertetalJ EoOOa). In the previous equation, we assume 
that the HpCO to H^CO ratio is equal to the isotopic ele- 
mental ratio, as sup ported by the available observations (e.g. 
ISchoier et alj f2002). Using this equation, the opacities values 
reported in Tableware obtained. 

The opacity values reported in Table @] are relatively low, 
which indicates that the lines are moderately thick, vwith the 
exception of the line at 351 GHz towards NGC1333-IRAS2. 
The uncertainty on the latter opacity is however relatively large, 
as shown by the errors bars reported in Table@] 
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Fig. 1. Spectra of the eigth H 2 CO transitions in TableEJobserved towards NGC1333-IRAS4A, NGC1333-IRAS4B, NGC1333- 
IRAS2, and L1448-MM respectively. 



4.2. Rotational diagram analysis 



method, and Man gum & Woottenlll993l for its application to 
formaldehyde lines). 

To obtain a first order estimate, we derived the beam-averaged Fig.0]shows the H2CO rotational diagrams of the observed 
column density of formaldehyde and rotational temperature sources. In these diagrams, the ortho to para ratio of formalde- 
by means of the standard rotational diagram technique (see hyde was kept as a free parameter, and was derived by min- 
iGoldsmifh & Langen Il999l for a general description of the imizing the x 2 between the observed fluxes and the rotational 
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L1448-N L1157-MM L1537 VLA1623 




15 -5 5 10 



J 1 1 Ml II 1 


1 1 1 1 1 1 _ 












V lar (km s ') 



V 1-r (km s ') 



V 1-r (km s ) 



5 5 10 



V,, r (km s ') 



Fig. 2. As Fig.EJfor the sources L1448-N, Ll 157-MM, L1527 and VLA1623. 



diagram predictions. The best agreement is obtained for a value 
of about 2 on all the sources. The fact that this value is lower 
than the high-temperature limit of 3 suggest that the formalde- 
hyde is formed at low temperature, around 20 K jKahane et alJ 
1984). However, while this ratio seems to be the same for all 
sources, we emphasize that it is highly uncertain. In particu- 
lar, the fact that the same transitions have been observed on all 



the sources can lead to systematic errors on this value. A more 
accurate derivation of the ortho to para ratio would need a cor- 
rection for the line opacities, which has the effect of scattering 
the p oints in the rotational diagram (see iGoldsmith & La nger 
I1999L for a review on the effect of opacities in a rotational 
diagram). This correction is not possible here as only a lim- 
ited number of H^CO transitions has been observed. The 52,4- 
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Table 2. Integrated fluxes of H2CO lines in T m b dV units . Upper limits are given as 2 cr. The "-" symbol indicates lack of the 
relevant observation. 







o-H 2 CO 








p-H 2 CO 






Transition 


2l,2-ll,l 


4l.4-3l,3 


5l.5-4l.4 


3o.3-2o,2 


32,2-22,1 


5o,5-4o,4 


5 2 ,4-42,3 


54,2 -44,1 11 


E u (K) 


21.9 


45.6 


62.4 


21.0 


68.1 


52.2 


96.7 


234 


A„j (s-») 


5.4 x If)" 5 


6.0 x 10- 4 


1.2 x 10- 3 


2.9 x 10- 4 


1.6 x 10- 4 


1.4 x 10- 3 


1.2 x 10- 3 


5.0 x 10- 4 


y (GHz) 


140.839 


281.527 


351.769 


218.222 


218.476 


362.736 


363.946 


364.103 


gu 


5 


9 


11 


7 


7 


11 


11 


11 


Telescope 


IRAM 


IRAM 


ICMT 


IRAM 


IRAM 


ICMT 


JCMT 


ICMT 


HPBW (") 


17 


9 


14 


11 


11 


14 


14 


14 


?7 m b or 6 c ff/ F c ff 


0.74 


0.47 


0.63 


0.62 


0.62 


0.63 


0.63 


0.63 


NGC1333-IRAS4A 


9.1 ± 1.4 


10.6 ± 2.6 


5.5 ± 1.2 


9.3 ± 1.9 


2.2 ± 0.4 


2.9 ± 1.0 


1.2 ±0.6 


1.7 ± 0.9 


NGC1333-IRAS4B 


6.8 ± 1.0 


12.1 ±3.0 


7.5 ± 1.7 


9.6 ± 1.9 


4.7 ± 1.0 


5.9 ±2.1 


3.6 ± 1.0 


0.9 ± 0.6 


NGC1333-IRAS2 


4.3 ± 0.6 


5.8 ± 1.5 


1.6 ±0.4 


4.9 ± 1.0 


1.0 ±0.2 


1.8 ±0.6 


0.6 ± 0.4 


<0.4 


L1448-MM 


3.3 ± 0.7 


4.7 ± 1.1 


1.0 ±0.2 


3.4 ± 0.6 


0.4 ±0.1 


1.3 ±0.4 


0.2 ±0.1 


<0.1 


L1448-N 


7.9 ± 0.9 


6.9 ± 1.6 


3.8 ±0.9 


5.7 ± 0.9 






0.6 ± 0.2 




L1157-MM 


1.2 ±0.2 


2.1 ±0.5 


1.2 ±0.3 


1.1 ± 0.3 


<0.2 


0.5 ± 0.2 


< 0.1 


<0.1 


L1527 


2.8 ± 0.7 


4.5 ± 1.1 


1.0 ±0.3 


3.0 ± 1.3 


0.2 ±0.1 


0.4 ± 0.2 






VLA1623 


3.8 ± 1.2 


2.7 ± 1.2 


0.9 ± 0.2 


5.0 ± 1.5 




1.2 ±0.4 


< 0.2 


<0.3 



Blended with the 5 4 ,i-4 4 _o H 2 CO line 



NGC1333-IRAS4A NGC1333-IRAS4B NGC1333-IRAS2 L1448-MM 
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Fig. 3. Spectra of the observed Hi 3 CO lines. 



42,3 line was only detected towards NGC1333-IRAS4A and to keep the derived parameters comparable from one source to 
NGC1333-IRAS4B. Because of its high energy (234 K), this the other. 

line is probably excited in hotter regions, and would increase Table summarizes the derived total column densities and 
the derived rotational temperature. This line was not included rotational temperatures. The column densities range from 2 to 

7 X 10 13 cm 2 , and the rotational temperatures from 1 1 to 40 K. 
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Table 3. Integrated fluxes of H2 13 CO lines in T m b dV units. Upper limits are given as 2 <x. The "-" symbol indicates lack of the 
relevant observation. 







o-H 2 J CO 




p-H 2 3 CO 


Transition 


2i,2-li,i 


4l,4-3l,3 


5l.5-4l, 4 


52.4-4 2 .3 


E u (K) 


21.7 


44.8 


61.3 


98.5 


A uJ (s-1) 


1.5 x 10- 4 


1.7 x 10~ 3 


3.4 x 10- 3 


1.1 x IO- 3 


v (GHz) 


137.450 


274.762 


343.325 


354.899 


HPBW (") 


17 


9 


14 


14 


Telescope 


IRAM 


IRAM 


ICMT 


ICMT 


rj mb or S c ff/F cff 


0.74 


0.47 


0.63 


0.63 


NGC1333-IRAS4A 


0.4 ± 0.1 


0.3 ±0.1 


0.3 ± 0.2 


<0.1 


NGC1333-IRAS4B 






0.3 ±0.1 


<0.1 


NGC1333-IRAS2 


0.2 + 0.1 




0.4 ± 0.2 


<0.1 


L1448-MM 






<0.1 


<0.2 


L1448-N 


0.2 ±0.1 


<0.1 






L1157-MM 






< 0.1 


<0.1 



Table 4. H2CO lines opacities derived from the H2 CO observations. 







o-H 2 CO 




p-H 2 CO 


Transition 


2 1.2- 1 1,1 


41.4-3J.3 


5l,5-4l,4 


52,4-42.3 


NGC1333-IRAS4A 
NGC1333-IRAS4B 


1 +(r/ 


5 +1 1 


1 2 +91J 
09+13 


< 2 
<0.5 


NGC1333-IRAS2 
L1448-MM 


1 l + ° 




I'M 

5.5_ 2 2 
< 2 


<4 


L1448-N 


O4+11 

-0.3 


<0.1 






L1157-MM 




< 1.4 





The values are both only lower limits to the actual gas tempera- 
ture and column density, as the gas temperatures can actually be 
significantly higher in the case of non-LTE conditions, and the 
derived column density can also be higher in case of optically 
thick emission. To correct for this effect, the derived column 
densities were recalculated adopting the average value of opac- 
ities quoted in Table@] The corrected column densities are also 
reported in Table|5| and range from 0.8 to 2 xlO 14 cm -2 . 



4.3. LVG modeling 

In order to derive the physical conditions of the emitting gas 
under non-LTE conditions, the formaldehyde emission has 
been modeled usin g an LVG code 5 . The collisional coefficients 
from lGreeiJlll99lh and the Einstein coefficients from the JPL 
database ( Pic kett et all 1 19981) were used. The LVG code has 
three free parameters: the column density to line width ratio 
N< ^&f 0) (which regulates the line opacity), the gas temperature 
r gas , and the molecular hydrogen density n(H2). When the lines 
are optically thin the line ratios only depend on the latter two 
parameters. Since we measured only marginally optically thick 
lines, the gas temperature and density were first constrained 
based on the line ratios predicted in the case of optically thin 
lines. The absolute line fluxes predicted by the model were 



5 D etails on the used LVG code can be found in ICeccarelli et all 
(2002). 



then compared with observations to constrain the H2CO col- 
umn density. 

The gas temperature and density have been determined by 
minimizing the^ 2 ed , defined as: 



^red 



1 



N-2 



(Observations - Model) 



(6) 



where all the observed H2CO lines were included. Fig.[5]show 
the x 1 contours. The derived r gas and n(H2) are reported in 
Table|5] The H2CO column densities were then constrained us- 
ing the observed 0-H2CO 5 i,5-4x 4 line flux, under the assump- 
tion of optically thin emission. 

The derived temperatures are between 30 and 90 K. These 
temperatures are significantly higher than the rotational tem- 
perature, suggesting that the observed transitions are subther- 
mally populated. Indeed, the derived densities support this con- 
clusion, as they vary between 1 and 6 x 10 s cm -3 . It is re- 
markable that the density range is so small, but the density es- 
timate is certainly biased because of the choice of the transi- 
tions. In particular, the densities are slight ly lower than those 
found by brake et all \ 1995b and Ivan Dishoeck et alJ dl995h 
for NGC1333-IRAS4A, NGC1333-IRAS4B, and IRAS16293- 
2422. This is probably due to the fact that those studies in- 
cluded only the higher frequency lines whereas we here in- 
cluded also lower frequency, and thus lower critical density 
lines. This re-inforces the conclusion that a range of densi- 
ties are present in the envelope, as predicted by the power- 
law density structure derived from continuum observations 
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50 100 150 50 100 150 50 100 150 

E„ p (K) E up (K) E up (K) 




Fig. 4. H2CO rotational diagrams derived for the observed sources. Circles and squares mark the ortho and para H2CO transitions 
respectively. Fluxes of the para transitions have been multiplied by the ortho to para ratio, obtained by minimizing the^- 2 between 
the observations and the predictions of the rotational diagram (see text). Solid lines show the best fit curves. 



Jj0rgensen et all2002l) . Finally, the H2CO column densities de- 
rived with the rotational diagram method are typically lower by 
less than a factor 5 than the ones derived using the LVG method 
with the exception of NGC1333-IRAS2. 

Table [5] also reports estimates of the average H2CO abun- 
dance in each source, obtained dividing the H?CO b y the 
H2 column densities derived bv I.Te>rgensen et all 1 2£)02) from 
submillimeter continuum observations. The latter refer to the 
amount of material with a temperature larger than 10 K, typi- 
cally at a distance of several thousands of AUs, more than the 
envelope encompassed by the beam of our observations. This 
material is likely an upper limit to the column density of the 
gas emitting the H2CO lines, but it provides a first approxi- 
mate estimate of the H2CO abundance. Table |5] shows rather 
large variations in the H2CO abundance from source to source. 
Given the approximation of the method used to derive them, 
this spread may not be entirely real. In the next section we an- 



alyze the observed lines by means of an accurate model that 
takes into account the source structure. 



5. Protostellar envelope model 

5.1. Model description 

The model used computes the line emission from a spherical 
envelope. Dust and gas have density and temperature gradi- 
ents, that are approximated as follows. The density profile is 
described by a power law n(Hi) °c r~ a , where a is between 0.5 
and 2. The case a = 1.5 corresponds to an entirely free-falling 
envelope, whereas a = 2 would mimic an isothermal sphere 
in hydrostatic equilibrium. The densities and dust temperature 

profiles of all the sources of the sample have been derived by 

i r 1 1 I r J 

J0rgensen et al. (2002), except for L1448-N, whose analysis is 
reported in AppendixlAl 
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Table 5. Results of the rotational diagram and LVG analysis. 







Rotational Dia 


*ram 






LVG 








Source 


Trot 


AWfHjCO) 

iniiiv z s 


A'CHjCO) 3 


r gas 


n(H 2 ) 


Af(H 2 CO) b 


H 2 CO/H 2 c 


CO/H 2 


1 




(K) 


(cm -2 ) 


(cm" 2 ) 


(K) 


(cm 4 ) 


(cm -2 ) 








NGC1333-IRAS4A 


25 


7 x 10 13 


2 x 10 14 


50 


3 x 10 s 


1 x 10 13 


5 x 10 -10 


8 x 10" 




NGC1333-IRAS4B 


40 


7 x 10 13 


2x 10 14 


80 


3 x 10 5 


2x 10 14 


7 x lO -10 


1 x 10" 


-5 


NGC1333-IRAS2 


24 


3 x 10 13 


1 x 10 14 


70 


3 x 10 5 


5 x 10 13 


1 x 10-'° 


2 x 10- 


-5 


L1448-MM 


19 


2 x 10 13 




30 


3 x 10 5 


6 x 10 13 


4 x 10- 10 


4 x 10- 


-5 


L1448-N 


22 


5 x 10 13 


8 x 10 13 


90 


1 x 10 5 


3 x 10 13 


7 x lO -10 






L1157-MM 


18 


1 x 10 13 




80 


3 x 10 5 


4 x 10 13 


3 x lO -11 


6x 10- 


-6 


L1527 


16 


3 x 10 13 




30 


6 x 10 5 


4x 10 13 


1 x 10- 9 


4 x 10- 


-5 


VLA1623 


11 


3 x 10 13 




80 


1 x 10 5 


8x 10 13 


3 x lO -10 


2x 10- 


-4 



a Corrected for opacity effects, assuming a value of t = 1 for NGC1333-IRAS4A, NGC1333-IRAS4B and NGC1333-IRAS2 and r = 0.4 

for L1448-N respectively (see Table|4] 
b Averaged over a 10" beam. 

c From H? column densities q uoted bv lj0rgensen et al] i2002h 
d From J0 raensen et alj i2002h 





Fig. 5.x 2 contours as a function of the density and temperature of the emitting gas. Contours indicates the 1,2 and 3cr confidence 
levels respectively. 
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The gas temperature profi le has then been computed by 
using the model developed by ICeccarelli et alJ Jl996l) . which 
solves the thermal balance in the envelope. In order to compute 
the gas temperature, one needs to solve the radiative transfer of 
the main coolants of the gas, i.e. H2O, CO and O. For the water 
abundance in the inner and outer regions we used the values 
derived by the analysis of IR AS 16293-2422 JCeccarelli et alJ 
l2000al) and NGC1333-IRAS4 dMaret et al.l2002l) : 4 x 1(T 6 and 
4 x 10~ 7 respectively. The CO abundance in the outer region 
is taken to be 10~ 5 (e.g. Llcirgensen et aTll2Q02l) . lower than the 
canonical value, as this species is depleted in the cold parts of 
the envelope. Finally, the oxygen abundance is taken to be 2.5 x 
10~ 4 . Fig.[6]shows as a typical example the case of NGC1333- 
IRAS2. The gas temperature tracks closely, but not completely, 
the dust temperature. In the very inner region the gas is colder 
than the dust because of the increase of the water abundance 
caused by the icy mantle evaporation (when Td ust > 100 K), 
which increases the gas cooling capacity. In the very outer re- 
gion the gas is colder than the dust too, because of the effi- 
cient gas cooling by the CO lines (see als o the discussion in 
ICeccarelli et all 19 96 and M aret et al.l2002h . These differences 
concern however small regions in the envelope, and therefore 
the results would be essentially the same if the gas and dust 
temperature are assumed to be equal. To fully quantify this ef- 
fect, we ran a model for NGC1333-IRAS2 with a gas temper- 
ature equal to the dust one, and no significant differences were 
found. 

The velocity field, which regulates the line opacity in the 
inner envelope, has been approximated as due to free falling 
gas towards a 0.5 M central object in all sources (no turbulent 
broadening is taken into account). In view of the importance of 
ice evaporation, the formaldehyde abundance across the enve- 
lope has been modeled by a step function: when the dust tem- 
perature is lower than the ice mantle evaporation (T dust ^ 100 
K) the abundance has the value X out , whereas it increases to 
X m in the 7d us t > 100 K region. Finally, a H2CO ortho to para 
ratio of 3 was assumed 6 . We will discuss the dependence of 
the obtained results on these "hidden" parameters in the next 
section. 

Finally, the line emission is computed by solving the radia- 
tive transfer in presence of warm dust in the escape probability 
formalism where the escape probability [3 is computed at each 
point of the envelo pe by integrating the f ollowing function over 
the solid angle Q. dCeccarelh et all 1994 : 



P = 



kL + kd (k L + k A ) 2 



- exp [- (k h + k d ) ALa] 
AL fh 



(7) 



where k^ and k d are the line and dust absorption coefficients 
respectively, and AL tn is the line trapping region, given by the 
following expressions: 



AL, 



th 



2Avth I - 



1 - -lC 

2 H 



(8) 



6 Since the value of 2 obtained with the rotational diagrams is 
highly uncertain, we prefer to adopt the canonical value of 3 for the 
standard model. We will show in section l5'.3.3l that the results are not 
substantially affected by this value. 



10 6 r 




100 



r (AU) 

Fig. 6. Density (solid line), gas (dashed line) and dust (dotted 
line) temperature across the envelope of NGC1333-IRAS2. 



in the infalling region of the envelope (where arcos (fj) is the 
angle with the radial outward direction) and 



AU =r 1- 



(9) 



in the static region (where R eiw is the envelope radius). In the 
previous equations, Avth is the thermal velocity width and v 
is the infall velocity. In practice, the photons emitted by the 
dust can be absorbed by the gas and can pump the levels of the 
formaldehyde molecules. At the same time, photons emitted 
by the gas can be absorbed by the dust (at the submillimeter 
wavelengths the dust absorption is however negligible). 

5.2. Results 

In order to constrain the inner and outer formaldehyde abun- 
dance in the envelope, we run a grid of models, varying X out 
between 10~ 12 and 10~ 8 , and X m between 10~ 12 and 10~ 4 re- 
spectively for each source. Both parameters were constrained 
by minimizing the x^ ed - The best fit model for each source was 
then obtained from the x^ ed diagrams shown in Fig.0 and the 
parameters are summarized in Table [S] The list of predicted 
0-H2CO spectra for each source are reported in AppendixlEl 

X out is well constrained in all sources, and varies between 
8 x 10~ n and 8 x 10~ 10 . These values differ on average by a 
factor 3 from the abundances derived by the LVG analysis. In 
four sources (NGC1333-IRAS4B, NGC1333-IRAS2, L1448- 
MM and L1527) the value of X^ a is also well constrained by 
the observations, with a 3cr confidence level. In three sources 
(NGC1333-IRAS4A, L1448-N and L1157-MM) we also de- 
tected formaldehyde abundance jumps, but the level of confi- 
dence is less or equal to 2cr. VLA1623 is the only source where 
no evidence of a jump is detected, although the data would not 
be inconsistent with it. The measured X m values are between 
1 X 10" 8 and 6x 10~ 6 , and the jumps in the formaldehyde abun- 
dance are between 100 and 10 4 . 
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Fig. 7. Protostellar envelope model /j^ ed contours as a function of the outer and inner H2CO abundances. The contours levels show 
the 1, 2, and 3 <x confidence levels respectively. 

Table 6. Formaldehyde abundances as derived from the envelope model in the outer parts of the envelope (X out ) and the inner 
parts of the envelope (X m ). The table also includes the radius where the dust temperature is 100 K and 50 K, and the density at 
the radius where the dust temperature is 100 K. 



Source 


^100K 


^50K 


n 


00K 


X 


GUI 






X 


n 








(AU) 


(AU) 


(cm 


- 3 ) 


















NGC1333-IRAS4A 


53 


127 


2 


X 


10 lJ 


2 


X 


10- 


10 


2 


X 


10- 


s 


NGC1333-IRAS4B 


27 


101 


2 


X 


10 s 


5 


X 


10- 


10 


3 


X 


10- 


6 


NGC1333-IRAS2 


47 


153 


3 


X 


10 s 


3 


X 


10- 


10 


2 


X 


10- 


-7 


L1448-MM 


20 


89 


2 


X 


10 s 


7 


X 


10- 


10 


6 


X 


10- 


-7 


L1448-N 


20 


95 


1 


X 


10 s 


3 


X 


10- 


10 


1 


X 


10- 


6 


L1157-MM 


40 


105 


8 


X 


to 8 


8 


X 


10- 


1 1 


1 


X 


10- 


-8 


L1527 


20 


140 


3 


X 


10 6 


3 


X 


10- 


10 


6 


X 


10- 


6 


VLA1623 


13 


48 


2 


X 


10 8 


8 


X 


10- 


10 










IRAS 16293-2422" 


133 


266 


1 


X 


10 s 


1 


X 


10- 


4 


1 


X 


10- 


7 



From lCeccarelli et alJ |2000c ) . 
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Fig. 8. Ratio of the model predictions over the observed fluxes 
of H2CO lines as a function of the upper energy of the line, 
for NGC1333-IRAS2. In the left panel a jump in the abun- 
dance (Tab.|6jl is adopted, while the in the right panel a constant 
abundance across the envelope is assumed. The circles and the 
squares represent H^CO and H^CO lines respectively. The ar- 
rows represent lower limits. 

To illustrate the reliability of the derived jumps, Fig. [8] 
shows the ratio between the model and the observations in the 
cases of no abundance jump and with a jump, for NGC1333- 
IRAS2 as an example. The constant abundance model can re- 
produce the observed fluxes of the lower transitions well, but 
it definitively underestimates the flux of the lines originating 
from the higher levels. 

In the next section, we discuss the effects of the other model 
parameters on the H2CO abundance determination. 

5.3. Dependence on other parameters of the model 

The derived formaldehyde abundances depend on four hid- 
den parameters: the adopted velocity and density profiles, the 
H2CO ortho to para ratio and the evaporation temperature. In 
the following we discuss the influence of all these parameters 
on the determination of the H2CO inner abundance. 

5.3.1 . Velocity profile 

In our analysis, we assumed a velocity profile of a free falling 
envelope, given by: 

(2GM\ 1/2 

v(r) = —J (10) 

where G is the gravitational constant and M the mass of 
the central star. The choice of a free falling velocity profile 
seems natural, as these protostars are believed to be in ac- 
cretion phase ("e.g. I Andre et alJfeoOOl) . Yet, the central mass is 
a p oorly constrained param eter. Recently, iMaret etall (|2002) 
and lCeccarelli et al . (2000a) derived a central mass of 0.5 and 
0.8 M Q for NGC1333-IRAS4A and IRAS 16293-2422 respec- 
tively. Here we adopt a central mass of 0.5 Mq for all the ob- 
served sources, but this parameter could vary from one source 
to the other. 

A different choice for the velocity profile could change the 
derived abundance. In particular, a higher central mass would 



L1448-MM .Y E red=3 I 

Static N = 10 




io-= 10- 7 io-= 10- 5 

Fig. 9. x 2 kA diagrams of L1448-MM for three different veloc- 
ity profiles. The upper left and upper right panels show the 
X 2 ed contours derived for a free falling envelope with a central 
mass of 0.3 and 0.8 M0 respectively. The lower panel shows 
the static envelope case, with a 1 km s _1 turbulent velocity. 
Contours indicate the 1, 2 and 3<x confidence levels respec- 
tively. 

imply a higher velocity gradient in the central parts of the enve- 
lope, and as a consequence, a lower opacity of the H2CO lines. 
This lower opacity would decrease the formaldehyde abun- 
dance needed to reproduce a given flux. These differences are 
expected to affect mainly the high lying lines, which originate 
in the inner parts of the envelope. 

In order to quantify this effect on the derived formaldehyde 
abundances, two models of LI 448 -MM were run, using a cen- 
tral mass of 0.3 and 0.8 M0 respectively. A third model was 
also run to check the case of a static envelope, where a turbu- 
lent line broade ning of 1 km s~ l has b een added. This model 
was adopted by J0rgensen et al . (2002), and reproduced well 
the observed low J CO emission, which originates in the static 
envelope. 

Fi gure I5~3 . 1 1 shows the;^ 2 ed diagrams obtained in the three 
cases. While the derived outer formaldehyde abundance is not 
much affected by the adopted velocity field, the inner abun- 
dance changes significantly when adopting a static envelope 
rather than a free-fall profile. Yet, the inner abundance is well 
constrained in the first two cases (free-fall with different central 
masses), and very weakly depends on the adopted central mass: 
8 x 10~ 7 and 5 x 10~ 7 for 0.3 and 0.8 M© respectively. On the 
contrary, only a lower limit is obtained if a static envelope is 
adopted: > 3 x 10~ 8 . This is due to different line opacities in 
the three cases. In a static envelope, the high lying lines become 
more easily optically thick, because of the reduced linewidth 
with respect to the free-fall cases. For this reason, these lines 
do not depend on X^ , when X[ n is ~ 10~ 7 , because they be- 
come optically thick. This explains why only a lower limit on 
X m can be computed in that case. 
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5.3.2. Density profile 

The H2CO abundances depend on the density profile derived 
from the continuum observations. In particular, the observa- 
tions used to derive the physical structure of the envelopes 
have been obtained with a typical beamwidth of 10", i.e. 2200 
AU at the distance of NGC1333. The observed continuum is 
therefore not very sensitive to the physical conditions in the in- 
nermost regions of the envelope, at scales smaller than a few 
hundred AUs. Consequently, the derived density power-law in- 
dex reflects rather the density in the outer part of the envelope, 
and the inner density, extrapolated from these power law, may 
be a rough estimate. Finally, the determination of the density 
profiles of some of the sources of the sample was difficult be- 
cause of the existence of mul tiple components ( J0rg ensen et alJ 
20021. [torgensen et al.l J2002I) reported an average uncertainty 
of +0.2 on the power-law index. If the density at the outer 
radius of the envelope is assumed to be correctly determined 
by the continuum observations, the uncertainty on the power- 
law index corresponds to an uncertainty of a factor five on the 
density at the inner radius of the envelope. In order to check 
the effect of this uncertainty on the derived abundances, we 
ran models of NGC1333-IRAS4B with an inner density arti- 
ficially multiplied by a factor 5 (note that the outer density is 
not changed). Whereas the H2CO outer abundance remains un- 
changed, the inner abundance decreases by about the same fac- 
tor 5. Uncertainties in the inner density could therefore lead to 
uncertainties on the derived inner abundances of the same order 
of magnitude. 

5.3.3. The H 2 CO ortho to para ratio 

The derived formaldehyde abundances depend also on the 
value of the H2CO ortho to para ratio. Given the relatively 
low number of observed lines, this parameter cannot be con- 
strained by the present observations. Actually, it is very badly 
constrained even in the case of IRAS 16293 -2 422, where many 
more fo rmal dehyde lines have be en observed. ICeccarelli et alJ 
(2000c) and ISchoier et al ] J2002I) report a value for the ortho 
and para ratio around 3, but with a large uncert ainty. We thus 
adopt ed the canonical value of 3 in our analysis dKahane et alJ 
1984). 

As an example, Fig.^|shows the influence of this parame- 
ter on the derived H 2 CO abundance of NGC1333-IRAS4A. We 
ran models with the ortho to para ratio 1 , 2 and 3 respectively. 
While the derived inner and outer abundances are almost iden- 
tical for the three ratios, the abundance jump is slightly better 
constrained for a ratio of 1 (3cr) than an higher ratio (2<x). 

5.3.4. Evaporation temperature 

Finally, the evaporation temperature r ev , at which the 
formaldehyde desorption occurs, is also a hidden parame- 
ter of our model. As described in §5.1, in the present study 
we adopted r ev — 100 K, which c orresp onds to the evapo- 
ration of water ices ( Aika wa et alJll99 7P. However, part of 
the desorption could also occur at the evaporation temper- 
ature of pure formaldehyde ices (50 K; Aika wa et al.l ll997: 



NGC1ooo-IRAS4A ;/ a red =2.7 
o/p = 1 N=12 




NGC1333-IRA34A X 3 red = 2.0 

o/p = 2 K=12 




10-' 10- 9 io- B 10- 7 10- 10 io-» io-= 10- 7 



KGC1333-IRAS4A x a red =2.3 
o/p = 3 N=12 




10-10 10-8 lo~« 10-' 



Fig. 10. xi d diagrams of NGC1333-IRAS4 for three different 
ortho to para ratios. The upper left, upper right and lower pan- 
els shows the xLa contours derived for an ortho to para ratio 
of 1, 2 and 3 respectively. Contours indicates the 1, 2 and 3cr 
confidence levels respectively. 



NGC1333-IRAS2 x™d=l-l 
T^=100K N=ll 




NGC1333-IRAS2 X r e d=2 3 




Fig. 11. x, ed diagrams of NGC1333-IRAS2 for an evapora- 
tion temperature of 100 K (left panel) and 50 K (right panel). 
Contours indicates the 1, 2 and 3cr confidence levels respec- 
tively. 



Rodge rs & CharnlevlEool) . For example, a detailed analysis 
of the formaldehyde spatial emission in IRAS 16293-2422 has 
shown that the formaldehyde abundance has a first jump, of 
about a factor 10, where T^st ^ 50 K , and a second jump o f 
about a factor 25 where r^st > 100 dCeccarelli et alJEoOlh . 
However, given the relatively small number of observed lines 
and the absence of spatial information on the formaldehyde 
emission in the source sample of the present study, we limited 
the H2CO abundance profile to a single step function. In order 
to check if the data are also consistent with a jump at 50 K, 
a model with a jump in the abundance at 50 K for NGC1333- 
IRAS2 was run (see Fig. II li . While X out is very little sensitive 
to this parameter, X m is about ten times smaller when assuming 
a jump at 50 K (2 x 10~ 7 and 2 x 10~ 8 for 100 and 50 K re- 
spectively). We note, however, that the best agreement with the 
data is obtained for an evaporation temperature of 100 K (xL* 
= 1.1 against 2.3 respectively). 
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5.3.5. Summary 

In the xL. analysis of §5.2, we adopted the most reasonable 
values of the various hidden parameters in the model. In this 
section we have explored the effect of variations in them on 
the derived abundance jumps. We emphasize that the actual 
value of the jump as well as the precise location are quite uncer- 
tain. Based on the previous analysis, the adopted velocity field 
seems to affect most strongly the amplitude and/or the presence 
of the derived jump, in particular when a static, turbulent field 
is considered. The other three parameters, the inner density, 
the ortho to para ratio, and the evaporation temperature seem 
to play a minor role. In this sense, the model assumption of an 
infalling envelope is probably the most critical in the present 
analysis. As already mentioned, we favor the hypothesis of col- 
lapsing envelopes, both because evidences are accumulating in 
this direction (e.g. lDi Francesco et al boo lb . and because it is 
the most natural one. 

To summarize, the sources (i.e. NGC1333-IRAS4B, 
NGC1333-IRAS2, L1448-MM and L1527) where the X 2 Kd 
analysis yields 3<r evidence for jumps, the presence of an 
abundance jump is rather firm in our opinion. Although more 
marginal, the data are consistent with the presence of a jump 
in the other surveyed sources as well. Appendix 151 lists the 
predicted fluxes of the brightest ortho formaldehyde lines. 
Predictions of para H?CO line fluxes c an be found on the web 
site MEPEW 7 dCeccarelli et al.ll2003l . In particular, the sub- 
millimeter lines are sensitive to the presence and amplitude of 
the jump in the H2CO abundance, and future observations with 
existing (JCMT, CSO) and future (e.g. SMA, ALMA) submil- 
limiter telescopes will better constrain this value. 

6. Discussion 

The first remarkable and by far the most important result of this 
study is the evidence for a region of formaldehyde evaporation 
in seven out of eight observed sources. In these regions, the 
formaldehyde abundance jumps to values two or more orders of 
magnitude larger than the abundance in the cold outer envelope. 
The transition is consistent with the location where the dust 
temperature reaches 100 K. The radius of these warm regions 
is between 13 and 133 AU, and the densities 8 are between 1 and 
20 x 10 s cm 4 . A straightforward interpretation of this result 
is that the grain mantles sublimate at 100 K, releasing into the 
gas phase their components, and, among them, formaldehyde. 
In addition, recent observations have shown th e presence of 
comp lex molecules, typical of massive hot cores dCazaux et alJ 
2003) towards IRAS 16293-2422, the first studied hot core of 
low m ass protostars llCeccarelli et alJl2000albld : ISchoier et alJ 
120021) . The similarity with the well studied hot cores of the 
massive protostars is certainly striking: warm, dense, and com- 
pact regions chemically dominated by the mantle evaporation. 
Even though the chemistry can be, and very probably is differ- 

7 http : //www- laog . obs . uj f- grenoble . f r/~ceccarel/mepei 

8 L1527 is an exception, with a density of 3x 10 6 cirr 3 , but, as com- 
mented bv lj0rgensen et all 120021) this may be due to the contribution 
of the disk, that may dom inate the conti nuum emission in the inner 
parts of this source lLoinard et al. 2002b). 
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Fig. 12. Derived H2CO X m abundances as a function of the den- 
sity (left panel) and the radius (right panel) where Td us t =100 
K. The star represent IRAS 16293-2422. 
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Fig. 13. Derived H2CO X m abundances as a function of 
ismm/iboi- The star represent IRAS 16293-2422. 

ent in high and low mass protostars, the hot cores represent ba- 
sically the same process in both type of sources. Our study does 
not address the possibility that some of the H2CO ice mantles 
may be liberated by shocks interacting with the inner envelope, 
since the line wings have been excluded from our analysis. The 
role of shocks could be studied by future high angular resolu- 
tion maps of the line center and line wing emission. 

6.1. X- m versus L smm /L bol 

The H2CO abundance in the inner region shows a loose trend 
with a variety of source characteristics. These include the den- 
sity in the inner region and the size of the region (Fig. I12I . as 
well as the ratio of the submillimeter luminosity to the bolo- 
metric luminosity, L smm jL\, \ (Fig. ll3t . In assessing these cor- 
relations, we should, of course, keep the large uncertainties - an 
order of magnitude - as well as systematic effects in mind. In 
particular, underestimating the density will immediately result 
r /rfte]Jt^ r htfe atm § tne abundance required to explain the obser- 
vations (cf. section l5~3l . Likewise, if we underestimate the size 
of the emitting region, we will need a higher H2CO abundance. 
Now, which of these, if any, correlations is the prime correla- 
tion and which one is derived is not directly clear. In particular, 
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within the inside-out-collapse, the density and inner radius are 
anti-correlated as are the luminosity and the radius of the evap- 
orated region. 

The apparent anti-correlation between the H2CO abun- 
dance and L smm ILf, i (Fig. 1 1 31 is of particular interest. The 
latter is generally taken as an indicator of the evolutionary 
stage of the protostar where a larger value for L smm ILb i im- 
plies "youth" (e.g. large amounts of cold material surround- 
ing the YSO). The anti-correlation may seem a surprising re- 
sult, as the most accepted scenario predicts that formaldehyde 
is formed on th e grain surfaces, likely by successive hydro- 
genation of CO jTielens & Hagenlll982tlcharnlev et aDl 19971) 
during the pre-stellar phase. Once in the gas phase because 
of the evaporation of the grain mantles, formaldehyde is de- 
stroyed (i.e. converted into more complex molecules) on a 
timescale of ~ 10 4 yr JCharnlev et aljfl99^) . In this picture, 
the youngest sources should have the largest X{ n , which is evi- 
dently not the case. This picture, however, might be somewhat 
over-simplified. Indeed, the process of ice evaporation is con- 
tinuous, involving larger and larger regions as the time passes 
and the luminosity of th e protostar inc reas es - as pointed out by 
the mo dels bv lCeccarelli et alJ ( fl9 96) and Rodge rs & Charnlevl 
J2003I) - so that the result is a shell of continuously refurbished 
H2CO, moving outwards. The main point is that the region of 
formaldehyde overabundance never disappears, but just moves, 
and the jump in the H2CO abundance is only governed by the 
quantity of formaldehyde in the grain mantles. 

If the L smm /L\, \ ratio is not an age indicator but rather a 
parameter affected more by the different initial conditions of 
each protostar, and specifically it is larg er in sources wher e 
the pre-stellar density is larger ("e.g. ljavawardhana et all200ll) . 
the trend of Fig. ^] would imply that the H2CO ice abun- 
dance depends on the pre-stellar conditions. Less dense pre- 
stellar conditions (i.e. lower L smm /Lboi ratios) would give rise 
to more H2CO enriched ices, because there is more atomic 
H and thus more grain surface hydrogenation to form H2CO. 
This is indeed consistent with the fact that the efficiency of 
CO hydrogenation into H2CO on the grain mantles is ~ 250 
times larger in HoO-ric h ices when compared to CO-rich ices 
JCeccarelli et alJ KOOll) . And since less dense regions have 
likely more H20-rich than CO-rich ices, because CO-rich ices 
likely form in relatively dense condensations (Bacmann et al. 
2002), the larger H2CO abundance in sources with a lower 
L smm /Lboi would therefore be due to a real larger efficiency of 
the H20-rich ices in forming H2CO. As a consequence, our 
finding would suggest that the L smm /Lboi ratio does not probe 
the evolutionary status of protostars, but rather mainly reflects 
their different initial conditions. 

Alternative explanations are also possible. For example, re- 
cent laboratory works suggest that the formation of formalde- 
hyde by CO hydrogen ation on the grains depends on the dust 
temperature ( Watanabe et al. 2003), and this may also be con- 
sistent with "older" protostars (i.e. lower L smm /Lboi ratios), 
being also hotter, having larger H2CO abundances. This, of 
course, would imply that the bulk of the H2CO is formed in 
a stage later than the CO c ondensation, namely during the 
pre-stellar core phase (Bac mann et al.lE()02t) . Whether this is 
likely is difficult to say, for CO may indeed be trapped on 



the grain mantles and partly converted into H2CO only when 
the grain temperature increases, as suggested by the labora- 
tory experiments. Another possibility is that formaldehyd e for- 
mation on grains needs U V radiation (e.g ld'Hendecourt et alJ 
1986; Schutte etal]ll996l) . Analogously to above, formalde- 
hyde would be formed only in a later stage, and the older the 
protostar, the larger the UV field and the larger the H2CO abun- 
dance. 

All these interpretations need to be taken with caution, of 
course, since the inferred variations in the inner H2CO abun- 
dance might reflect uncertainties in the density and/or size of 
the region emitting the H2CO lines. A similar study on a larger 
sample and focussing on higher energy lines is required to draw 
more definitive conclusions. 



6.2. Low versus high mass protostars 

In section 13 we have examined the evidence for the presence 
of jumps in the H2CO abundance in the warm gas surround- 
ing low mass YSOs. The presence of such abundance jumps 
in hot cores around high mass stars is not well established. On 
the one hand, the prototype of hot cores in regions of massive 
star formation - the hot core in the Or ion BN/KL region - has 
a H 2 CO abundance of 10~ 7 JSutton et al.lll995h . On the other 
hand, in a study of hot c ores in a sample of massive protostars, 
Ivan der Tak et alJ J2000I) did not find evidence for the presence 
of H2CO abundance jumps, but did find evidence for jumps 
in the CH3OH abundance. It is unclear at present whether the 
Orion hot core or the van der Tak sample is more representa- 
tive for the composition of hot cores in regions of massive star 
formation. Presuming that the differences in H2CO abundance 
jumps are real, we note that the composition of the ices - which 
drive the chemistry in hot cores - may well differ between re- 
gions of low mass and high mass star formation. 

Supporting this idea, the deuterium fractionation is dra- 
matically different in the high and low mass protostars. 
Doubly deuterated formaldehyde and methanol have been ob- 
served to be extremely abundant in low mass protostars when 
compared to massive protostars. The D2CO/H2CO ratio is 
more than 25 times larger in low than in high mass proto- 
stars JCeccarelli et alJll998l lLoinard et alJl2002al) . Deuterated 
methanol may be as abundant than its main isotopqmer i n 
the low mass protostar IRAS 16293-2422 JParise et alJ | 2002 ). 
whereas it is only 4% of CH3OH in Orion Jjaca et alJ ll993). 
Since thi s extreme deuteration is likely a grain mantle prod- 
uct (e.g. lCeccarelli et alJ200lHParise et all2002l) . the dramatic 
differences in the molecular deuteration between low and high 
mass protostars are already a very strong indication that man- 
tles in both type of sources are substantially different. This in- 
deed does not comes as a surprise, as the pre-collapse condi- 
tions very likely differ: warmer in high than in low mass stars, 
at the very least. 



6.3. X out versus CO abundance 

Finally, Fig. 1 141 compares the H?CO a bundance X om with the 
CO abundance derived bv LF0rgensen et alJ J2002B . in the outer 
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Fig. 14. Derived outer formaldehyde abundance X oul as a func- 
tion the CO abundance of the outer enveloppe. Circles rep- 
resent Class sources, The star represents IRAS 16293-2422, 
the square represe nts VLA1623, and the t riangles represent the 
pre-stellar cores of iBacmann eTaD J20021 1 2003 ) sample. 

envelope. On the same plot we also reported the values foun d 
in the prestellar cores studied bv lBacmann et alj ( 120021 12003). 
The first thing to notice is the similarity of the values found 
in Class sources and pre-stellar cores, in both molecules, 
despite the different methods used to derive the abundances. 
The similarity of the values adds support to the thesis that 
the pre-stellar cores of the Bacmann et al. sample are pre- 
cursors of Class sources, and that the outer regions of the 
envelopes of Class sources are formed by material still un- 
changed by the collapse, i.e. that reflects the pre-collapse con- 
ditions. Second aspect to note of Fig.^^is the clear correlation 
between the H2CO and CO abundance. In this case the interpre- 
tation is straightforward: in the outer, cold envelope molecules 
are depleted, and the degree of depletion is th e same for the 
CO an d the H2CO molecules. As discussed in Bac mann et alJ 
il2002l) . the limited CO depletion observed in pre-stellar cores 
strongly suggests that a relatively efficient mechanism (cos- 
mic rays?) re-injects CO into the gas phase. Since the bind- 
ing energies of the CO and H2CO are relatively different (e.g. 
Aika wa et alJll997h . the linear correlation of Fig. El strength- 
ens t he claim that H9CO molecules are trapped into CO-rich 
ices JCeccarelli e t al. 2001). 

7. Conclusions 

We have presented a spectral survey of the formaldehyde emis- 
sion of a sample of eight Class protostars, carried out with the 
IRAM 30m and JCMT telescopes. A total of eight transitions 
were selected for each source, covering a large range of upper 
level energies in order to probe different physical conditions. 
When possible, the isotopic lines were observed in order to de- 
rive the line opacities. Most of the lines are relatively narrow 
with a small contribution from wings extending to larger ve- 
locities. Using the standard rotational diagram method, we de- 



rived rotational temperatures between 1 1 and 40 K, and H2CO 
column densities between 1 x 10 13 and 7 x 10 13 cm -2 . For 
the sources with detected H^CO lines, opacity corrections in- 
crease the derived column densities to the range 0.8 and 2 
xlO 14 crrT 2 . In order to test the effect of non-LTE excitation, 
the observed emission has been modeled using a LVG code. 
In this way the derived temperatures are significantly higher 
than the rotational temperatures, suggesting that the observed 
transitions are sub-thermally populated. The inferred densities, 
between 1 and 6 x 10 5 crrT 3 , indeed support this conclusion. 

To take into account the density and temperature gradients 
in the protostellar envelopes, the emission has been modeled 
using densities and dust temperature profiles derived from pre- 
vious studies of the continuum emission of these objects. The 
gas temperature in the envelopes was computed using a code of 
the thermal balance for protostellar envelopes. The formalde- 
hyde abundance across the envelope has been approximated by 
a step function : an outer abundance X out where 7d us t ^ 100 K, 
and a inner abundance X[ n at 7d us t S 100 K. We show that the 
outer abundance, X oul , is well constrained in all the sources, 
and varies between 8 x 10" 11 and 8 x 10~ 10 . The inner abun- 
dance X m is well constrained in four sources with a 3 cr level 
confidence, while in three sources it is only a < 2 <x evidence, 
and no evidence of a jump is found in VLA1623. The derived 
values of X m range between 1 x 10~ 8 and 6 x 10~ 6 , showing 
jumps of the formaldehyde abundance between 2 and 4 orders 
of magnitude. The most important conclusion of this study is 
hence that large amounts of formaldehyde are injected into the 
gas phase when the grain mantles evaporate at 100 K. 

We have discussed the uncertainties on the actual values of 
the hidden parameters of the model, namely the velocity and 
density profile, the H2CO ortho to para ratio, and the evapora- 
tion temperature. The uncertainty in these parameters makes 
the abundance jump value and jump locations uncertain for 
some sources. Future observations of higher frequency lines 
and modeling of the line profiles may distinguish between the 
different interpretations. 

We found that sources with lower L smm /Lboi ratios pos- 
sess the largest inner H2CO abundances. We discussed why we 
think that this reflects very likely a different H2CO enrichment 
of the grain mantles. 

We found that the inner H2CO abundances are systemat- 
ically larger than the H2CO abundances of the hot cores of 
the sa mple of massive protostars studied by Ivan der Tak et all 
(2000). This supports to the idea that low and high mass proto- 
stars have a different grain mantle composition. 

Finally, the derived outer H2CO abundance shows a 
clear correlation with the CO abundance, implying that both 
molecules are depleted by a similar factor in the outer envelope, 
namely that H2CO molecules are likely trapped in CO-rich ices 
in the outer envelope. 
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Table A.l. L1448-N best fit parameters from the DUSTY mod- 
elling and derived physical parameters. 



DUSTY parameters 

y iToo 

a 1.2 
r IO o 0.8 



Derived physical parameters 





10 AU 


RlOK 


7000 AU 


VR, 


4.9 x 10 s cm" 3 


nioK 


1.1 x 10 5 cm" 3 


M, 0K 


2.5 M Q 



Appendix A: Density and temperature profile of 
L1448-N 

The density and temperature profile of L1448-N have been determined 
following the method used by J0rgens en et alj i2002). The spectral en- 
ergy distribution (SED) and JCMT-SCUBA maps at 450 and 800 ftm 
have been compared to the prediction of the radiative code DUSTY 9 
llvezic & ElitzuJl997l) . The JCMT observations were taken from the 
JCMT archive. Fig. IA. II shows the result of the fits of the brightness 
profile and SED of this source. The envelope parameters from this 
modeling are summarized in Table lA. fl 

L1448-N shows a relatively flat density profile. This profile may 
reflect the asymmetry of the source and t he exten sion of the emission, 
which can flatten the derived profile (see Jorge nsen et alj2002l for a 
discussion of the effects of asymmetries in the derived density profile). 

Appendix B: Models predictions for formaldehyde 
lines fluxes 

In this Appendix we present the model predictions for the fluxes of 
selected 0-H2CO transitions. In Table |B~T1 we report the line fluxes 
computed assuming a constant abundance X ollt across the envelope, 
while Table IB .2*1 reports the fluxes predicted assuming an abundance 
X oM in the outer part of the envelope where Tdust < 100 K, increasing to 
X in in the inner part of the envelopes, at the radius where T iust > 100A*. 
The adopted values of X m and X out are the best fit values reported in 
[5| In these two tables, the line fluxes are expressed in erg s~' cm 4 . 
Although a conversion in K km s -1 would have been more practical 
to compare these fluxes with observations, this conversion depends on 
the telescope used and the extent of the line emission. However if the 
emitting region is smaller than the telescope beam, the fluxes can be 
approximat ively converted into K k m s -1 using the conversion factors 
reported by Ceccarel li et alJ 1200 3l) . 

While low lying lines fluxes are comparable in two tables, higher 
energy lines fluxes substantially differ. Higher frequency observations 
can therefore help to distinguish between a constant abundance in the 
envelope and a jump in the abundance, especially on the sources of 
our sample for which the jump in the abundance is uncertain. 



9 DUSTY is publicly available at 

http : / /www . pa . uky . edu/~moshe/dusty/ 
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Fig. A.l. Observations and result of the modelling of L1448N. The upper left panel shows the SCUBA observations at 450 and 
850 yum. The upper right panel shows the SED (diamonds) and the fitted model (solid line). The lower left panel shows the 
brightness profile and the fitted model. The dashed line represents the beam profile. The lower right panel shows the derived 
temperature dust (solid line) and density profile (dashed line). 
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Table B.l. Predicted line fluxes of selected 0-H2CO transitions for a constant abundance X out across the enveloppe. 



Transition 


Freq. 
(GHz) 


(CITT 1 ) 








Fluxes 
(erg s -1 cirr 2 ) 








IRAS4A 


IRAS4B 


IRAS 2 


T 1/1/1QA/TA/T 
bl44olvlJVl 


T 1/1 A CM 
L14451N 


L1157MM 


LI 527 


VLA1623 


2 1,2 - 


1 1,1 


140.8 


15 


2.6E-15 


1.6E-15 


4.6E-16 


1.1E-15 


2. IE- 15 


3.7E-16 


7.6E-16 


3.8E-16 


2 u _ 


1 

1 1.0 


150.4 


16 


2.5E-15 


1.4E-15 


3.9E-16 


n oc 1 A 
9.zb-lo 


11C K 

2. lb- 15 


3.5E-16 


6.4E-16 


3.6E-16 


3l,3 _ 


2 1,2 


211.2 


22 


3.4E-15 


1.6E-15 


6.3E-16 


1 1 c 1 c 
l.lb-lj 


*2 1.C 1 £ 


5.6E-16 


5.7E-16 


4.3E-16 


3 1,2 - 


2 1,1 


225.6 


23 


2.8E-15 


LIE- 15 


4.6E-16 


A ^c 1 A 
O.Ob-lo 


O AC 1 ^ 

z.4b-0 


4.4E-16 


2.6E-16 


3.2E-16 


4l,4 - 


3l,3 


281.5 


32 


3.3E-15 


LIE- 15 


6.6E-16 


H ?C 1 A 


O 1 C K 

2. lb-lj 


5.3E-16 


2.0E-16 


3.2E-16 


43,2 - 


3 3,1 


291.3 


98 


3.0E-17 


3.3E-18 


9.7E-18 


A AC 1 O 

4.0b-lo 


OC 1 
Z.ob-15 


3.7E-18 


3.1E-19 


1.5E-18 


43,1 - 


^3,0 


291.3 


98 


3.0E-17 


3.3E-18 


9.7E-18 


A AC 1 O 

4.0b-l5 


O OC 1 

2. ob- to 


3.7E-18 


2.9E-19 


1.5E-18 


4l,3 - 


Q 

J 1,2 


300.8 


33 


2.6E-15 


5.8E-16 


4.8E-16 


A 1 C 1 A 

4.lb-lo 


1 1 c 1 c 

i.ib-ij 


3.9E-16 


5.6E-17 


2. IE- 16 


5 i,5 - 


4 1,4 


351.7 


43 


2.7E-15 


5.3E-16 


6.0E-16 


A AC 1 A 

4.4b- 10 


OIC 1 A 

o.ib-10 


4.0E-16 


4.9E-17 


2.0E-16 


C 

~>3,3 _ 


/I 

43,2 


364.2 


110 


8.6E-17 


7.1E-18 


2.6E-17 


^c 1 Q 


j.lb-lo 


1.0E-17 


2.9E-19 


3.7E-18 


C 

■5 3,2 — 


/I 

43,2 


364.2 


110 


8.6E-17 


7.1E-18 


2.6E-17 


r» cc 1 


£ i c 10 
3.ib-i» 


1.0E-17 


2.7E-19 


3.7E-18 


c 

3 1.4 _ 


/I 

4 1,3 


375.8 


46 


2.5E-15 


3.8E-16 


5.2E-16 


nc 1 A 
j.Ub-lo 


c f\c 1 A 

j.Ub-10 


3.6E-16 


1.6E-17 


1.6E-16 


Ol,6 - 


c 

J 1,5 


421.9 


57 


2.3E-15 


3.1E-16 


5.6E-16 


o cic 1 A 
z.9b-lo 


Q /1C 1 A 

3.4b- to 


3.2E-16 


1.4E-17 


1.4E-16 


z: 

03,4 - 


C 

->3,3 


437.1 


125 


1.6E-16 


1.0E-17 


4.7E-17 


1 AC H 

1.4b-l / 


A /1C 1 O 

o.4b-lo 


1.8E-17 


2.2E-19 


5.9E-18 


03,3 - 


c 

J3,2 


437.2 


125 


1.6E-16 


1.0E-17 


4.7E-17 


1 ac \n 
1.4b-l / 


A /1C 1 O 

o.4b-lo 


1.8E-17 


2.0E-19 


6.0E-18 


A 

"1,5 - 


c 

3 1.4 


450.8 


61 


2.3E-15 


2.3E-16 


5.1E-16 


01c 1 a 
z.lb-lo 


O OC 1 A 

2.2b- to 


3.0E-16 


3.9E-18 


1.2E-16 


7i,7 - 


<c 

□1,6 


491.9 


74 


2.0E-15 


1.9E-16 


5.1E-16 


1 nc 1 A 

i.yb-10 


1 AC 1 A 

1.0b- 10 


2.7E-16 


3.9E-18 


1.0E-16 


*7 

'3,5 ~ 


f. 

03,4 


510.1 


142 


2.5E-16 


1.2E-17 


6.8E-17 


1 nc 1 n 
1. /b-1 / 


A AC 1 O 

o.ob- 18 


2.7E-17 


1.4E-19 


7.7E-18 


'3.4 - 


03,3 


510.2 


142 


2.5E-16 


1.2E-17 


6.8E-17 


1 nc 1 n 
1. /b-1 / 


A AC 1 O 

o.ob- 18 


2.7E-17 


1.3E-19 


7.7E-18 


<n 

'1,6 - 


"1,5 


525.6 


78 


2.0E-15 


1.5E-16 


4.9E-16 


1 AC 1 A 

Lob- lo 


1 1 C 1 A 

l.lb-lo 


2.6E-16 


1.2E-18 


9.0E-17 


Q 

°1,8 ~ 


*7 

'1,7 


561.8 


93 


1.8E-15 


1.2E-16 


4.7E-16 


1 AC 1 A 

1.4b-lo 


nc 1 n 
o.Ub-1 / 


2.2E-16 


1.2E-18 


7.4E-17 


Q 

03,6 - 


*7 

'3,5 


583.1 


161 


3.4E-16 


1.3E-17 


9.1E-17 


1 OC 1 1 

l.yb-1 / 


A OC 1 O 

o.zb-18 


3.6E-17 


8.2E-20 


8.9E-18 


O 

03,5 _ 


'3.4 


583.2 


161 


3.4E-16 


1.3E-17 


9.0E-17 


1 CC 1 T 

1.5b-l / 


A OC 1 O 

o.zb-18 


3.6E-17 


7.6E-20 


8.8E-18 


O 

81,7 " 


'1,6 


600.3 


68 


1.9E-15 


1.1E-16 


4.7E-16 


1 OC 1 A 

l.zb-lo 


A 1 C 1 "7 

0. ib-1 / 


2.3E-16 


4.3E-19 


6.9E-17 


5*1,9 - 




Ol,8 


631.6 


114 


1.7E-15 


8.7E-17 


4.4E-16 


1 1 C 1 A 

l.lb-lo 


/i nc 1 -7 

4. /b-l / 


2.0E-16 


4.4E-19 


5.8E-17 


n 

9 3 7 _ 


O 

03,6 


656.1 


183 


4.4E-16 


1.4E-17 


1.1E-16 


1 CIC 1 O" 


C CC 1 O 


4.5E-17 


4.4E-20 


9.5E-18 


^3,6 


83,5 


656.4 


183 


4.4E-16 


1.4E-17 


1.1E-16 


1 OC 1 T 

l.yb-1 / 


c cc 1 


4.4E-17 


4.1E-20 


9.5E-18 


9i, 8 - 


8l,7 


674.7 


121 


8.0E-17 


1.8E-15 


4.6E-16 


9.8E-17 


3.9E-17 


2. IE- 16 


1.8E-19 


5.6E-17 


lOi.io 


-9 l9 


701.3 


137 


6.5E-17 


1.7E-15 


4.2E-16 


8.3E-17 


2.9E-17 


1.9E-16 


1.7E-19 


4.6E-17 


10 3 ,8 - 


-9 3 , 7 


729.1 


207 


1.3E-17 


5.3E-16 


1.3E-16 


1.8E-17 


4.3E-18 


5.1E-17 


1.9E-20 


9.2E-18 


10,,, - 


-9i, 8 


729.6 


207 


1.3E-17 


5.2E-16 


1.3E-16 


1.8E-17 


4.3E-18 


5.1E-17 


1.8E-20 


9.2E-18 


lll.ll 


- 10 U o 


749.0 


146 


6.2E-17 


1.8E-15 


4.4E-16 


7.9E-17 


2.5E-17 


1.9E-16 


7.7E-20 


4.5E-17 


lluo 


- 10i,io 


770.8 


163 


5.0E-17 


1.7E-15 


4.0E-16 


6.5E-17 


1.9E-17 


1.7E-16 


6.7E-20 


3.7E-17 


H1.10 


- 10l, 9 


823.0 


173 


4.9E-17 


1.8E-15 


4.3E-16 


6.5E-17 


1.7E-17 


1.9E-16 


3.4E-20 


3.7E-17 


12l,12 


- lll.ll 


840.2 


191 


3.6E-17 


1.6E-15 


3.7E-16 


4.8E-17 


1.1E-17 


1.6E-16 


2.2E-20 


2.8E-17 


12!,!! 


- liuo 


896.7 


203 


3.6E-17 


1.8E-15 


3.9E-16 


4.8E-17 


1.0E-17 


1.7E-16 


1.3E-20 


2.8E-17 
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Table B.2. Predicted line fluxes of selected 0-H2CO transitions for an abundance X out in the outer enveloppe, increasing to X m at 
the radius where 7d USt > 100 K. 



Transition 


Freq. 
(GHz) 


E up 

(cnr 1 ) 








Fluxes 
(erg s - ' citT 2 ) 








IRAS4A 


IRAS4B 


IRAS 2 


L1448MM 


L1448N 


L1157MM 


LI 527 


2 1,2 - 


1 1,1 


140.8 


15 


2.7E-15 


1.6E-15 


4.9E-16 


1.1E-15 


2. IE- 15 


3.7E-16 


7.6E-16 


2i,i - 


1 1,0 


150.4 


16 


2.5E-15 


1.4E-15 


4.2E-16 


9.2E-16 


2.1E-15 


3.6E-16 


6.4E-16 


3l,3 - 


2 1,2 


211.2 


22 


3.5E-15 


1.7E-15 


7.3E-16 


1.1E-15 


3.3E-15 


5.9E-16 


5.7E-16 


3 1,2 - 


2i,i 


225.6 


23 


2.9E-15 


1.1E-15 


5.8E-16 


6.9E-16 


2.4E-15 


4.8E-16 


2.7E-16 


4l,4 - 


3l,3 


281.5 


32 


3.5E-15 


1.1E-15 


8.9E-16 


8.1E-16 


2.2E-15 


6.2E-16 


2.1E-16 


43,2 - 


3 3,1 


291.3 


98 


1.6E-16 


9.7E-17 


2.8E-16 


6.4E-17 


4.6E-17 


4.5E-17 


9.2E-19 


43,1 - 


33,0 


291.3 


98 


1.6E-16 


9.6E-17 


2.8E-16 


6.4E-17 


4.6E-17 


4.5E-17 


9.0E-19 


4i, 3 - 


3l,2 


300.8 


33 


2.9E-15 


6.8E-16 


7.7E-16 


4.9E-16 


1.1E-15 


5.1E-16 


6.0E-17 


5 1.5 - 


A 

4 1,4 


351.7 


43 


3.1E-15 


6.9E-16 


1.0E-15 


5.6E-16 


9.2E-16 


6.2E-16 


5.7E-17 


J3.3 - 


A 

43,2 


364.2 


110 


4.6E-16 


1.9E-16 


5.6E-16 


1.5E-16 


1.1E-16 


1.5E-16 


2.3E-18 


-> 3,2 — 


A 

43,2 


364.2 


110 


4.6E-16 


1.9E-16 


5.6E-16 


1.5E-16 


1.2E-16 


1.5E-16 


2.3E-18 


5l,4 - 


4 1,3 


375.8 


46 


3.0E-15 


5.7E-16 


1.1E-15 


4.5E-16 


6.4E-16 


6.2E-16 


2.5E-17 


61,6 - 


5 1,5 


421.9 


57 


3.0E-15 


5.7E-16 


1.3E-15 


4.9E-16 


5.4E-16 


7.2E-16 


2.7E-17 


°3,4 - 


c 

~>3,3 


437.1 


125 


9.2E-16 


3.2E-16 


9.6E-16 


2.6E-16 


2.1E-16 


3.2E-16 


3.9E-18 


63,3 - 


53,2 


437.2 


125 


9.2E-16 


3.2E-16 


9.6E-16 


2.6E-16 


2.2E-16 


3.2E-16 


4.0E-18 


"1,5 - 


c 

J 1.4 


450.8 


61 


3.1E-15 


5.5E-16 


1.4E-15 


4.7E-16 


4.6E-16 


7.9E-16 


1.7E-17 


7 1,7 " 


s: 

Ol,6 


491.9 


74 


3.1E-15 


6.0E-16 


1.7E-15 


5.2E-16 


4.6E-16 


9.2E-16 


2.0E-17 


7 3 ,5 - 


63,4 


510.1 


142 


1.5E-15 


4.9E-16 


1.5E-15 


3.9E-16 


3.4E-16 


5.7E-16 


4.8E-18 


'3,4 - 


03,3 


510.2 


142 


1.5E-15 


4.9E-16 


1.4E-15 


3.9E-16 


3.4E-16 


5.6E-16 


4.9E-18 


7 1,6 - 


61,5 


525.6 


78 


3.4E-15 


6.6E-16 


1.9E-15 


5.6E-16 


4.8E-16 


LIE- 15 


1.7E-17 


81,8 _ 


7 1,7 


561.8 


93 


3.4E-15 


7.5E-16 


2.2E-15 


6.3E-16 


5.4E-16 


1.2E-15 


1.8E-17 


83,6 


7 3,5 


583.1 


161 


2.3E-15 


7.4E-16 


2.2E-15 


5.9E-16 


5.2E-16 


9.1E-16 


4.9E-18 


8 - 


7 3,4 


583.2 


161 


2.3E-15 


7.4E-16 


2.2E-15 


5.9E-16 


5.2E-16 


9.1E-16 


5.0E-18 




01,7 - 


7 1,6 


600.3 


68 


3.9E-15 


8.6E-16 


2.6E-15 


7.2E-16 


6.2E-16 


1.5E-15 


1.5E-17 


9l,9 - 


Q 

°1,8 


631.6 


114 


4.1E-15 


1.0E-15 


3.0E-15 


8.3E-16 


7.2E-16 


1.7E-15 


1.5E-17 


"3,7 


03,6 


656.1 


183 


3.2E-15 


1.0E-15 


3.1E-15 


8.3E-16 


7.3E-16 


1.3E-15 


3.9E-18 


9 3 ,6- 


83,5 


656.4 


183 


3.2E-15 


1.0E-15 


3.1E-15 


8.3E-16 


7.3E-16 


1.3E-15 


4.0E-18 


9i, 8 - 


81,7 


674.7 


121 


4.7E-15 


1.2E-15 


3.6E-15 


9.6E-16 


8.5E-16 


2.0E-15 


1.3E-17 


lOi.io 


-9,,, 


701.3 


137 


5.1E-15 


1.4E-15 


4.1E-15 


1.1E-15 


9.8E-16 


2.3E-15 


1.0E-17 


10 3 ,8 - 


-9 3 , 7 


729.1 


207 


4.2E-15 


1.4E-15 


4. IE- 15 


1.1E-15 


9.2E-16 


1.8E-15 


2.2E-18 


10,,, - 


-9i, 8 


729.6 


207 


4. IE- 15 


1.4E-15 


4. IE- 15 


1.1E-15 


9.1E-16 


1.8E-15 


2.2E-18 


lll.ll 


- lOi.io 


749.0 


146 


5.8E-15 


1.6E-15 


4.8E-15 


1.3E-15 


1.1E-15 


2.7E-15 


8.5E-18 


H1.10 


- lOi.io 


770.8 


163 


6.3E-15 


1.8E-15 


5.4E-15 


1.4E-15 


1.3E-15 


3.0E-15 


6.4E-18 


H1.10 


- 10l, 9 


823.0 


173 


7.3E-15 


2.1E-15 


6.3E-15 


1.7E-15 


1.5E-15 


3.5E-15 


5.1E-18 


12l,l2 


- 1 1 1,1 1 


840.2 


191 


7.9E-15 


2.4E-15 


7.0E-15 


1.8E-15 


1.6E-15 


3.7E-15 


2.9E-18 


12 U1 


- Hl.10 


896.7 


203 


9.1E-15 


2.8E-15 


8.2E-15 


2.1E-15 


1.8E-15 


4.3E-15 


2.4E-18 



